BINS = 256/16;
BINS = 20;

mct=[176 21 6 278 30 18 19 18 202 16 5 341 37 17 295 8 209 130 11 456 43 22 5 261 43 616 464 13 400 20 8 89 130 11 38 129 8 10 393 80 732 11 15 210 66 4 194 55 14 124 107 7 3 70 88 13 524 10 190 10 298 7 5 364 17 11 3 10 2 305 491 7 832 119 8 10 325 321 14 206 405 62 810 101 6 6 270 509 465 168 8 8 7 5 15 14 9 9 3 10 24 34 125 9 11 126 820 671 17 7 15 94 351 9 103 6 13 14 515 11 25 10 7 9 60 18 486 124 593 1079 10 261 176 575 214 387 4 13 466 5 9 7 5 973 95 110 9 3 4 359 114 5 10 55 20 91 99 59 8 11 48 7 107 61 11 11 14 169 86 6 15 5 148 6 658 134 389 490 117 9 61 27 101 11 103 55 7 7 7 12 51 8 399 8 5 39 3 374 225 225 434 125 10 88 6 4 12 14 7 110 7 9 497 57 9 10 11 19 368 90 30 79 33 103 30 12 9 152 5 110 245 761 132 122 5 136 8 13 20 16 53 136 577 854 748 3 3 6 46 4 14 13 357 229 11 525];
%mct =ones(256)*450
mct=[81 91 89 105 88 93 103 73 112 91 94 99 104 88 103 94 106 105 98 80 103 88 98 98 94 100 113 95 102 95 113 92 95 97 106 105 115 96 97 98 106 95 92 89 100 108 90 98 116 119 98 89 95 85 103 76 105 96 92 107 97 88 116 86 72 102 89 95 102 102 93 94 93 87 92 101 105 100 85 99 96 93 90 92 97 94 77 94 80 88 100 109 103 117 94 120 95 100 87 82 101 113 108 107 85 110 110 93 111 86 96 98 112 109 70 127 94 90 120 97 116 109 79 106 104 103 98 104 92 80 87 96 89 125 93 114 93 111 90 116 99 99 99 112 102 97 68 86 105 90 87 100 98 103 87 113 94 88 107 118 118 103 115 109 84 93 92 100 83 90 93 106 89 82 95 100 91 69 111 110 113 108 101 95 98 80 91 104 93 87 87 100 102 95 86 100 98 108 98 90 93 99 105 89 91 96 120 96 85 82 108 123 97 106 104 114 107 96 85 98 106 108 95 88 98 82 85 89 103 110 102 76 101 100 119 89 110 103 84 83 101 104 95 94 93 92 89 85 100 82 73 91 100 98 74 116];
mct =[4 1 3 4 2 5 1 4 2 6 1 4 2 3 5 2 1 0 3 3 7 2 4 3 5 2 1 5 0 7 4 8 1 1 1 6 3 3 5 4 5 2 4 6 4 4 6 1 3 4 2 3 4 1 3 3 3 1 5 3 4 3 5 5 4 1 2 1 4 2 5 4 3 1 4 1 6 4 5 2 5 3 6 3 2 4 4 6 4 7 7 5 5 5 2 4 4 6 0 1 2 3 3 1 1 4 5 3 3 1 6 3 2 4 8 2 5 3 5 3 1 7 3 3 4 5 3 2 5 6 1 1 4 3 3 5 2 4 3 3 4 3 4 2 4 3 3 3 6 5 5 5 6 0 2 3 2 1 4 1 5 3 3 0 4 2 2 3 3 5 3 9 5 2 2 9 7 3 1 1 1 1 3 2 4 3 4 1 3 4 3 4 3 3 2 11 4 1 2 5 4 3 2 2 4 7 4 3 2 6 3 3 1 2 1 0 7 3 4 5 4 6 5 5 1 5 3 1 3 2 3 7 2 3 5 5 5 3 3 3 4 4 5 2 0 4 1 3 5 8 0 2 3 1 1 2 ]

xh = 1:BINS;
ch = zeros(BINS,1);

for i=1:BINS
    e = (xh(i)./BINS)*(4);
    xh(i)= 10^(-1+e);
end



for i=1:256
    f = log(mct(i))/log(10)+1;
    if f<0
        f = 0;
    elseif f>4
        f = 4;
    end
    
    j=floor(f/(4/BINS))+1;
    ch(j) = ch(j)+1;
end

semilogx(xh,ch,'r');

ch